Synthetic Gauge Fields for the Vibrational Excitations of Trapped ions 
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The vibrations of a collection of ions in a microtrap array can be described in terms of tunneling phonons. 
We show that the vibrational couplings may be tailored by using a gradient of the trap frequencies, together 
with a periodic driving of the trapping potentials. These ingredients allow us to induce effective gauge fields on 
the vibrational excitations, such that phonons mimic the behavior of charged particles in a magnetic field. In 
particular, microtrap arrays are well suited to realize a quantum simulator of the famous Aharonov-Bohm effect, 
and observe the paradigmatic edge states typical from quantum-Hall samples and topological insulators. 
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Introduction - The ultimate goal of quantum simulation 
(QS) is to find experimental platforms where quantum many- 
body physics can be explored [1]. This challenge requires 
efficient methods to prepare quantum states, an exquisite con- 
trol of the interactions, and precise measurement techniques. 
Quantum-information technologies have an important appli- 
cation in this context, since they provide us with a powerful 
toolbox for the manipulation of quantum systems. In partic- 
ular, trapped ions [2] are an interesting candidate, their main 
advantage being an unrivaled efficiency in preparing and mea- 
suring quantum states at the single-particle level. Besides, the 
strong long-range Coulomb interactions make them suitable 
for the QS of a variety of collective phenomena, from quan- 
tum magnetism [3] to dissipative models [4]. So far, there 
have been experiments with up to 9 ions [5], and a big effort is 
being focused on scaling them up. A promising avenue are the 
so-called two-dimensional arrays of microtraps (2DAM) [6], 
which may open new routes towards the many-body regime. 
Unfortunately, this setup still faces some hurdles, such as the 
large distances between ions (d x ^40 /im) leading to weak 
spin- spin interactions. In order to realize QS schemes based 
on vibration-mediated interactions, it is fundamental to over- 
come these issues. Alternatively, focusing directly on the vi- 
brational modes [8] yields a significant speed-up with respect 
to decoherence rates. In fact, the transfer of vibrational exci- 
tations between two aligned traps was recently observed [9] . 

In this Letter, we show how to tailor the vibrational cou- 
plings in a 2DAM, such that this speed-up is exploited. This 
opens the possibility of building a QS of lattice bosons under 
synthetic gauge fields. We note that laser-based methods for 
neutral atoms might also lead to effective gauge fields [10]. 
Our proposal, however, relies on the different concept of 
photon-assisted tunneling [11, 12], and requires a gradient of 
the individual trapping frequencies, together with a periodic 
driving of the trapping potentials that can be achieved by an 
optical force. This Letter is structured as follows: i) We show 
that the amplitude and phase of the vibrational couplings be- 
tween ions can be tuned by inducing resonances that corre- 
spond to the absorption/emission of photons from a classical 
driving field (photon-assisted tunneling), ii) We extend this 



result to 2D, and show how it leads to the implementation of 
synthetic gauge fields, where phonons move like charged par- 
ticles in a lattice. Hi) We present an implementation of the 
required drivings by means of optical forces, such that the 
optical phase can be intrepreted as an effective gauge field. 
iv) We propose a proof-of-principle of our ideas with four 
ions in a plaquette displaying a discrete version of the cele- 
brated Aharanov-Bohm effect [13]. v) We suggest to concate- 
nate those plaquettes in ladders, leading to Aharonov-Bohm 
cages [14], and allowing to observe the edge states character- 
istic of quantum-Hall samples and topological insulators [15]. 

i) Photon-assisted tunneling - We introduce our scheme for 
two ions with mass M and charge e, trapped by independent 
potentials with frequencies CO\ 2 (Figs. l(a)-(b)). The equi- 
librium positions, separated by d x , lie along the x direction, 
and the axial vibrational modes are periodically driven. The 
Hamiltonian is H (t) = H (z) +H C , with (ft = 1) 

e 2 

H (T) = Y,co j a)a j +H d (T), H c = -3 (5*i - 8x 2 f. (1) 

7=1,2 d x 

H c is the Coulomb coupling to second order in the ion dis- 
placements, 8xj = {aj + y/2M(Dj 9 with Oj(aj) phonon 
creation (annihilation) operators. The periodic driving is 

#d(0= £ «dCOs(co d T + 7 -)flta 7 -, (2) 
7=1,2 

where (0^ (rf^co^) is the driving frequency (strength), and (j>j a 
site-dependent phase. We assume that C0\ = CO, o>2 — G> + Aft), 
and {Aco, r\&(0&\ <C CO, namely, both the frequency difference 
and driving strength are small perturbations to the trapping 
frequency. In absence of driving, the vibrational coupling is 

H c =J c (a\a l + a[a 2 ), (3) 

where J c = -fico, and /3 = e 2 /Mco 2 d^. Equation (3) holds 
for |/ c | <C CO, such that the cross-terms, a\a 2 , a\a\, can be ne- 
glected in a rotating-wave approximation (r.w.a.). This con- 
dition, which is met for the experiments in [9], allows us to 
interpret the dynamics as the tunneling of phonons [8]. 
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Figure 1. Arrangement of ion micro traps: Schematic represen- 
tation of the microtrap layout for: a) two-ion link, c) four-ion pla- 
quette; e) many-ion rhombic ladder (the relevant parameters of the 
effective Hamiltonian (13) are also shown). Requirements for the 
photon-assisted tunneling of phonons: b) The frequencies of two ad- 
jacent traps ft>2 are shifted by Aco. By shining a pair of Raman 
lasers, one assists the phonon transfer; d) For a plaquette, the gradient 
is along x and the Ak has a component along the x-y-z axes. 



To understand the effects of driving, we express Eq.(3) in 
the interaction picture with respect to Hq(t), where 0/(t) = 



-/^•Tg-iTjdsiiiCfl^T+^JgiTjdsin^). After the trivial transforma- 



tion aft 



,t*T] d sin(0i) 



—> au one writes H c (t) by replacing the bare 
Coulomb coupling J c by a time-dependent dressed coupling 



,«((B d T+lte) -lV(O) d T+0i) 



(4) 



where ^f s ( r ld) are Bessel functions of the first kind. By 
choosing the driving frequencies r(Q& = Aco, with r = 1,2, 
one selects resonant processes that correspond to the absorp- 
tion/emission of r photons from the classical driving field. If 
Aco ^> / c , only the resonant terms must be considered after a 
r.w.a., and the effective vibrational coupling can be written as 



i§(0i+fc) 



^(T] d ,A0) =^^(T] d )^ +r (T] d )^ + 2 



)A0 



(5) 



where A(j) = — fa . Since no perturbative assumption is re- 
quired on the driving strength, the dressed coupling may 
be close to the bare one J c . Fig. 2a) shows a calculation of the 
dressed coupling under different conditions. By tuning r/d and 
{(j>j}, one controls the amplitude and phase of the tunneling, 
which may be enhanced or even completely suppressed. We 
note that the coherent control of tunneling is interesting on its 
own [11], and can now be investigated with trapped ions. 



ii) Synthetic magnetic fields- We extend this scheme to dif- 
ferent geometries given by the ion equilibrium positions in a 
2DAM, separated by d Xl d y , and labeled by vectors of integers, 
i = iy)- In the most general situation, trapping frequencies, 
co a i, depend on the site i and the spatial direction a = x,y,z. 
The trap potentials, together with a driving term, are 



(6) 



The vibrational couplings between ions in the array arise due 
to the Coulomb interaction 



1 



<5rr 



Srjl 



(7) 



where r?, and (<5rj) a = (a a ^ + a} a .) / y/2MG) a ^, are the equi- 
librium positions and relative ion displacements. We assume 
that the vibrational modes in different directions are not cou- 
pled, and the phonon number is conserved. The validity of 
this approximation is quantified below. In the harmonic ap- 
proximation, that is, up to second order in 5rj, we find 



i>j,a 



3W-j)«W-j) 



(8) 



IM^JCDaiCDaj 



|r? 



where r?_j = r? — r? . The assumption of independent vibra- 
tions in each direction holds for \co a ,i — &>/3,j|a^/3 ^ 
whereas phonon number conservation is valid if co a { ^> |/". . | . 

In quantum mechanics, charged particles under electromag- 
netic fields acquire a phase that depends on the field back- 
ground. To make a QS of this phenomenon, we focus on the 
ion motion in direction a, and choose a linear gradient along 
x > ®a,i = ®a + Ao)/ x , together with phases that depend linearly 
on the position, fa = (f) x i x + (jtyiy. Equation (2) is generalized 
to H&{i) = LiT]dCO ( iCos(ft) ( iT + (jh) a lc m analogy t° the 
two-ion case, we find that the effective vibrational couplings 
to leading order in 7]^ [16] are the following: 



+ 1 ' 



J c;i,j u ijcJjc' 



(9) 



where Afaj = fa — fa, and Si x j x is the Kronecker delta. The 
first term describes the photon-assisted tunneling along x, 
whereas the second one is the bare coupling along y. A crucial 
result is that the amplitude of tunneling around a plaquette, 



w a rot jcl rot rot 

YV J [r] [r] M M ;i+*v 



5 = |Wg|e^, 



yields an accumulated phase that depends on the laser param- 
eters 0o = —r(j)y, and can be recast in terms of the celebrated 
Aharonov-Bohm phase [13], <fc = e§ cbt • A, where § is the 
line integral along the plaquette, and A = r(j) y y / \ed x d y )x is 
a synthetic vector potential. Accordingly, phonons move as 
charged particles subjected to a magnetic field perpendicular 
to the microtrap array, and yield a bosonic counterpart of the 
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Azbel-Harper-Hof stacker model [17]. We stress that arbitrary 
fluxes 0o £ [0,2k) can be attained, even reaching one flux 
quantum per unit cell; a regime unaccessible in solid-state ma- 
terials for realistic magnetic fields. This opens the possibility 
to observe a dipolar version of the fractal Hofstadter butterfly, 
among other interesting effects presented in iv)-v). 

Hi) Realization of the periodic driving- The simplest setup 
to realize Eq. (2) would consist of an array of microtraps, 
where the driving fields are provided by the local control of 
the electrodes. Since this scheme is yet to be realized and 
scaled [6, 7], we base our alternative approach on state-of-the- 
art optical forces. We focus on the vibrational modes trans- 
verse to the microtrap plane, a = z, although other schemes 
along the x-y plane are equally valid [18]. We consider lasers 
that drive two-photon stimulated Raman transitions between 
the electronic levels of the ions |0)i, |l)i. The lasers are de- 
tuned by and provide a Raman wavevector Ak, 

Hl= ~t% 0i (e iAk(r i +5ri) "^ T + H.c.) , (10) 

where £2l is the Rabi frequency, and 0\ is an operator acting 
on the electronic levels. By a proper choice of the laser de- 
tunings and polarizations, one may realize 0\ = li, or other 
operators like 0{ = a? that widen the applicability of our QS 
(Outlook). The effect of the ion-laser interaction can be un- 
derstood as a periodic driving of the microtrap frequencies 
under the assumptions below. We consider a gradient along 
x, co z [ = CQ Z + Acoi x , such that the following conditions are 
fulfilled {Aco, <C {o) a , \co a - G)p\ a ^p}. Let us define the 
Lamb-Dicke parameter along a, r\ a = \Ak a \/y/2M(O a . In the 
limit 7] a <C 1, we perform a Taylor expansion of (10) up to 
second order in = //l,o +#l,2- Note that #l,o 

does not affect the vibrational modes, and can be ne- 
glected if QlTIoc <C co a in a r.w.a. This leads to 

#l,2-^l £ riaripO { cos(Ak-rf - co L z)al [ i a^ i . (11) 
ct,j3,i 

Finally, by considering \oo a — (dp\a^p ^ ^lVcc we ne- 
glect the coupling between different directions, and get the an- 
nounced periodic driving presented in Eq. (2) with the follow- 
ing identifications, co^ = r\&CD& = ^L^f , and <ft = — Ak r?. 

Current microtrap design [6, 7] is consistent with the above 
requirements, \co a - CQp \ a ^} > {cDl, Aco} >/y. Typ- 
ically, co a /27l « 1-10 MHz, J/271 « 5 kHz. To fit the inequal- 
ity, we can take 6)^/271 « 50 kHz. With a typical Lamb-Dicke 
parameter of r\ a « 0.2, the condition v\a^L <^ \co a — C0p | is 
still fulfilled [18]. In Fig. 2b), we compare the effective de- 
scription (5) to the exact optical forces (10) for a two-ion ar- 
ray, with parameters Aco = 0.05co z , r\ z = 0.2, £2 L = 0.75co z , 
j3 = 0.002, r = 1, where the phonon Hilbert space is truncated 
to rc max = 4. We observe an excellent agreement between both 
descriptions, yielding assisted tunneling for A(j) « 71 [11]. 

/vj-vj Aharonov-Bohm physics in lattices - We apply our 
ideas to a square lattice, and set r = 1, |Ak x |<i x = 2nn x + X , 



|Ak-y|t/y = 2^n-y + 0_y , with %, ^ e Z. The latter are intro- 
duced because typical ion distances are larger than optical 
wavelengths. From Eq. (9), we get the tight-binding model, 

#e ff = L/f 1];M+ Ww~''^+ I 4uW^i + »* + H - C -' 
i i,m>0 

(12) 

where <fo = (f) y [19]. Note that photon-assisted tunneling 
along the diagonals has been neglected, since for (j) x = 27l — (j) y 
that tunneling amplitude vanishes &{y\&,2it) = (Fig. 2a)). 
Besides, the remaining diagonal terms, i+jt+m ^, are negli- 
gible for m > 1 due to the fast dipolar decay. 

iv) Discrete Aharonov-Bohm effect - The simplest realiza- 
tion of this tight-biding model consists of a single plaquette 
(Figs. lc)-d)). In Figs. 2c)-d), we test the validity of the effec- 
tive dynamics in (12), by comparing with an exact numerical 
calculation of the complete driven Hamiltonian (10). These 
results describe a realization of the discrete Aharonov-Bohm 
effect with minimal required resources. In Fig. 2d), we ob- 
serve that an initial excitation can follow two possible paths, 
either 1^2^3orl^4^3, enclosing a net flux <fe = n. 
The paths interfere destructively and forbid the phonon to tun- 
nel to site 3. Conversely, in Fig. 2c), phonons tunnel around 
the plaquette for 0^ = 0. 

v) Aharonov-Bohm cages and flatband physics- Let us con- 
sider an interesting route beyond the single plaquette, which 
is the rhombic 3-leg ladder presented in Fig. Id). This system 
is described by the Hamiltonian 

tf=£/i(feja y + cj^ (13) 

j 

where we have labelled the boson operators for each leg as 
aj,bj,Cj (Fig. le)). This Hamiltonian follows directly from 
Eq. (12), when the plaquettes are arranged along a diago- 
nal, with Ji e 2 /(2mdy),J 2 h^r(j[&, A(j))(d y /d x ) 3 , and 
(j) = 0-y, X = 271 — (j) y . This model yields two effects. Due to 
the Aharonov-Bohm interference for (j) = ±7t, all the modes of 
the system are localized and one obtains flat vibrational bands 
(Fig. 2e)). In particular, the non-zero energy modes corre- 
spond to the so-called Aharonov-Bohm cages, where phonons 
are not allowed to tunnel two plaquettes apart [14]. Besides, 
one finds the so-called edge states, which are mid-gap modes 
exponentially localized around the boundaries. By tuning (j), 
and J2/J1, one can explore a transition between two topologi- 
cally non-equivalent phases (Fig. 2d). 

Finally, let us consider the experimental requirements for 
the implementation of our ideas. The duration of a QS to ob- 
serve the effects of the synthetic gauge fields is of the order of 
l/J[ r ], being of the order of the bare couplings J c , which 
are in the range 1-2 kHz [9]. This can be increased to 5 KHz 
following the trap design [7], and even enhanced by orders 
of magnitude by further miniaturizing the electrode structure, 
and storing more than one ion per lattice site [9]. The main 
competing decoherence mechanism is heating of the motional 
modes [2]. Heating rates as low as 0.07 phonons /ms have 
been reported in cryogenic traps [9], in principle allowing to 



4 





(H :«*>' 
|H-<X>-), 


M 


□□□□□□□□□ 


X 


h\0>-)\ 

>-\-<o)\ 



10 m 20 



31 




Figure 2. Photon-assisted tunneling and Aharonov-Bohm cages: (a) Effective hopping amplitude, \^ r= \\, as a function of (T]d,A0). 
The values of the dashed line at r/^ « 0.6 are used in (b). (b) Photon-assisted hopping for a two-ion link of a single vibrational excitation 
| Xf/o) =a[\0) under the effective description (5) (red line), and the complete Hamiltonian (3) with the driving term (10) (yellow dots). We plot 
the vibrational population n\ transferred to site 2 after time t* = 7r/|7j r ](0.6, 0) |. Maximum phonon transfer occurs at A(j) = %. (c,d) Evolution 
of the phonon excitation in a rectangular plaquette. We plot the phonon populations under approximation (9) (Pi(t)), and under the exact 
Hamiltonian Eqs. (6)-(8) driven by (10) (n,-(f)). In (c) (j) = 0,(j) x = n, (j) y = 0,d x = dy\^i(r} d , tt)| 2 / 3 , 

«max — 2, and other parameters same as 

(b). In (d), 0o = 71, and there is an Aharonov-Bohm destructive interference that inhibits tunneling to site 3. (j) x = 7i,(j)y = 71, £1^ = 0.25ft) z , 
and other values same as (c). (e) Energy spectrum E for the 7T-flux regime of Hamiltonian (13) for N = 31 microtraps. Flat bands appear at 
E « ±2J\ (Aharonov-Bohm cages) (see also the schematic description of the eigenstates), and also at E = (which arise solely due to the 
geometry of the lattice). Also, in the middle of the gaps, single edge states localized to the boundaries of the ladder arise, (f) Energy level 
spectrum as a function of the effective flux, which displays a gap-vanishing point at = 0. 



implement our ideas. Even when heating rates are compara- 
ble to couplings /[ r ] , they may induce a thermal background 
over which propagation of vibrational excitations may still be 
observed [9] . Note that experimental techniques are available 
for preparation and measurement of phonon states [20] . Also, 
the vibrational spectrum can be measured without local ad- 
dressing in the ions fluorescence sidebands [2]. 

Conclusions and Outlook - We have presented a proposal 
to induce synthetic gauge fields for ions in microtrap arrays, 
which is based on the photon-assisted tunneling of vibra- 
tional excitations. By considering trap designs with anhar- 
monicites, effective phonon-phonon interactions can be in- 
cluded [8], which may allow us to study strongly correlated 
phases. Inducing electronic state-dependent drivings, one gets 
effective spin-orbit couplings that induce disorder [21]. Also, 
by adding dissipation, i.e. motional heating, one may study 
quantum effects in energy transport in the presence of noise 
[22]. These ingredients make a versatile QS of many-body 
physics, which would outperform classical computers for « 10 
ions, and « 4 phonons per ion. That size seems feasible in the 
near future in view of current experimental progress [6]. Fi- 
nally, our scheme could be extended to other systems such as 
photons in arrays of cavities in circuit QED [23]. 
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